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Abstract 

In this paper we construct the Hamiltonian constraint operator of loop quantum cosmology using 
holonomies defined for arbitrary irreducible SU(2) representations labeled by spin J. We show that 
modifications to the effective semi-classical equations of motion arise both in the gravitational part 
of the constraint as well as matter terms. The modifications are important for phenomenological 
investigations of the cosmological imprints of loop quantum cosmology. We discuss the implications 
for the early universe evolution. 
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I. INTRODUCTION 



In the study of cosmology many questions remain unanswered despite the wealth of 
observational data at hand. In particular our understanding of the early universe is lacking. 
It is at these high energies that any classical cosmological models are no longer relevant and 
a quantum theory of gravity is needed. A deep quantum understanding of the regime is 
required to answer, for instance, questions related to the initial big-bang singularity and to 
understand what conditions preceded the onset of chaotic inflation. 

A particularly exciting quantum theory of cosmology is known as loop quantum cosmol- 
ogy (LQC) (for a review and references see The progress made in LQC to date has 
been impressive. LQC has provided explanations to long standing questions in quantum 
cosmology. The reason for its success can be traced to the fact that it has as its basis a 
candidate full theory of quantum gravity known as loop quantum gravity (for reviews see 
3 El 0, 5], (J). Thus LQC is expected to handle on a firmer basis the extreme high energy 
regime of the early universe where non-perturbative effects are paramount. For instance it 
has been shown that the theory contains a robust mechanism for singularity avoidance Q. 
The mechanism arises as a natural feature of the theory and does not rely on ad-hoc forms 
of matter. This is in contrast to earlier forms of quantum cosmology such as those based on 
Wheeler-DeWitt quantization where singularity avoidance does not arise naturally. 

Yet, an even more exciting result is that LQC predicts modifications to the evolution 
of the universe for scales larger than the Planck scale. For instance it has been shown 
that LQC can modify the behavior of the inverse volume such that it becomes a bounded 

n 

function This gives rise to modifications to the behavior of matter in the form of a 
scalar field which in turn predicts a period of super- inflation at energies near the Planck 
scale In this super-inflationary regime the inflaton is pushed up its potential due to the 
existence of an anti-friction term in the Klein-Gordon equation which can set the stage for 
the onset of slow-roll inflation . The modified dynamics in this regime can even lead 

to measurable signatures at large scales in the CMB 11]. 

Most of the phenomenological studies to date have exploited the modification of the 
inverse volume. For instance, in a collapsing universe with a scalar field the anti-frictional 
term in the Klein-Gordon equation becomes frictional and a bounce can occur Q]. This 
effect has been studied for a scalar field with a self-interaction potential in positive curvature 
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oscillatory universes where the effects of repeated expansion/contraction push the scalar field 
up its potential providing another mechanism to establish the initial conditions for slow-roll 
inflation Q| . In addition the modifications can lead to a bounce in the context of cyclic 
models with colliding branes, possibly removing singular behavior that has plagued these 
models 

These results are intriguing yet some open questions remain as to their validity. The 
modifications to the inverse volume arise by exploiting an ambiguity in defining the inverse 
volume operator which is required to quantize the matter part of the Hamiltonian constraint 
for a scalar field. In defining the inverse volume operator one traces over SU(2) valued 
holonomies and the ambiguity lies in choosing the irreducible representation in which to 
perform the trace. The ambiguity is labeled by a half-integer J representing the spin of the 
representation. A neglected fact is that this same ambiguity arises in the gravitational part 
of the Hamiltonian constraint where one also traces over SU(2) valued holonomies. So far, 
the gravitational part of the constraint has been quantized using the J = 1/2 representation 
while the ambiguity has been freely exploited for the matter part. 

It is therefore not clear if corrections arising in the gravitational part might change the 
dynamics and alter significantly the phenomenological consequences. From the full theory 
the gravitational part contains an inverse volume factor and the possibility remains that 
modifications here might cancel the effects in the matter part. It is thus of critical importance 
to determine the complete quantum corrections for arbitrary J arising in both terms of the 
Hamiltonian constraint. 

In this paper we aim to clear up the uncertainty with respect to this ambiguity. We 
explicitly construct the isotropic Hamiltonian constraint operator for arbitrary spin J rep- 
resentations. We show that the resulting difference equation is higher order for larger J and 
that geometric quantities in the gravitational part of the constraint obtain quantum modifi- 
cations analogously to the inverse volume modifications for the matter part. We show that 
the higher J constraint operator contains the same singularity resolution mechanism attest- 
ing to its robustness. We will derive and motivate effective semi-classical equations of motion 
from the constraint operator and comment on the phenomenological consequences and show 
that the qualitative results do not change. As an explicit example, the J = 1 constraint 
operator is constructed and analyzed with special attention payed to the semi-classical limit. 



3 



II. LOOP QUANTUM COSMOLOGY 



We start with a brief overview of LQC and its kinematical Hilbert space. In Dirac 
quantization for systems with constraints (such as general relativity) one first starts by 
constructing a kinematical Hilbert space by quantizing the system ignoring the constraints. 
The constraints are then represented as self-adjoint operators on this space and physical 
solutions are wave functions annihilated by the constraint operators. In this paper we 
will restrict ourselves to homogeneous and isotropic models with zero spatial curvature. 
Upon the imposition of this symmetry reduction, the Hamiltonian of isotropic GR consists 
of a single constraint, the Hamiltonian constraint. The details of the construction of the 



16]. 



We will follow the 



kinematical Hilbert space and the constraint operator are given in 
same conventions and notation here. 

The first step toward quantization is a classical symmetry reduction of the action of 
general relativity written in connection form. In the full theory of loop quantum gravity 
the action is written in terms of a connection A which contains the information about 
curvature, as well as a triad E conjugate to A which encodes the geometry. After reduction 
to homogeneity and isotropy the connection and triad are labeled by single quantities c and 
p respectively. These variables satisfy the commutation relation 

{C,P} = 

where k = 8tcG and 7 is the Barber-Immirzi parameter which represents a quantum ambi- 



guity (black hole entropy calculations can fix its value to 7 ~ .2735 







18| ) . In terms of 



the standard metric variables given by the scale factor a , the new variables are related as 

\p\ = a 2 (la) 
c = 7a (lb) 

Note the absence of the factor of 1/2 in the formula for the connection (llbj) which appears 
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The correct factor (without the 1/2) can be derived using the relation from the full 
theory A\ = T l a + 7 A"* where A % a is the full connection, K l a is the extrinsic curvature, and 
r* is the spin connections which vanishes identically for isotropic flat models. The absolute 
value around the triad component p indicates that we are extending the classical phase space 
of isotropic GR to include both orientations of the triad (positive and negative p). Note 
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that this does not mean that negative p corresponds to a negative scalar factor. For both 
regions the scale factor, and hence the volume, is positive. 

Dynamics is determined completely by the Hamiltonian constraint H given by 

H = — sgn(p) ^/\f\c 2 + H m (2) 

where H m is the matter Hamiltonian. Note again the factor of 3 which differs from previous 
works. With the correct factors here the resulting action Sgr = J dt^pc — NH GR is 
equivalent to the standard isotropic Einstein-Hilbert action Sqr — j^J y/\g\R, and the 
Hamiltonian equations p = {p, H} and c = {c, H} are equivalent to the equations of motion 
derived from the Einstein-Hilbert action. 

With the symmetry reduced action we now turn to quantization. The theory is quantized 
using ideas and techniques from the full theory of loop quantum gravity. The basic variables 
to be quantized are holonomies of the connection and the smeared triad integrated on a 
two-surface. Using this, the quantum configuration space inherits a discrete structure and 
an orthonormal basis is given by states labeled by a parameter \i given by 

(c\ fj) = exp(ific/2) fiGl. 

The basis states are eigenstates of the triad operator p with eigenvalues 

V |A*> = ^ M • (3) 

Geometric operators are built from the triad, which allows us to understand the physical 
meaning of the parameter \i. In particular it is an eigenstate of the volume operator V = 
\p\ 3 / 2 hence the physical volume is given by 

3/2 

(4) 
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More complicated operators are built from the basic triad operator p and holonomy 
operators. Because of homogeneity and isotropy we only need consider holonomies along 
straight lines which can be parameterized as 

hi = exp(fi cTi) 

where n are the generators of the Lie algebra of SU(2) satisfying [r,, tj] = e^r fc , and // is pro- 
portional to the length of the holonomy. In the fundamental representation the holonomies 



are given by the simple formula 

= cos(/i c/2) + 2Tj sin(yU c/2) 

where the notation (2) indicates that the holonomy is defined in the fundamental represen- 
tation (J = 1/2). In this paper we will use superscripts to the left of an object to indicate 
the irreducible representation of SU(2) used. The holonomy operators in turn behave as 
exponentiated position operators and thus act on the basis states by finite shifts. More 
precisely, the holonomy operators consist of the cosine and sine operators which we will for 
compactness notate cs and sn with actions 

cs = cos O^yJ I/ 1 ) = \ Ia* -»- A*o> |A* — Vto) 

sn \fi) ee sin ( ^- ] \fi) = -- + fi ) - |/i - /x ) (5) 



v 2 J " ' 2 

Knowing the action of the holonomy operators let us now turn to two essential operators: 
the inverse volume and the gravitational part of the Hamiltonian constraint. The inverse 
volume operator is needed to quantize certain forms of matter, for instance a scalar field. 
Since the eigenstates of the volume operator are normalizable, the naive inverse volume 
operator (with eigenvalues equal to the inverse of the volume eigenvalues) is not a densely 
designed self-adjoint operator. The solution is to recast the classical formula for the inverse 
volume using the Poisson bracket between the connection and the volume ^| . When quan- 
tized the connection is represented with a holonomy operator, and the Poisson bracket is 
replaced with a commutator. The resulting operator is diagonal in the basis and is given 
by H 
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^V- 1 \fi) 



yl/2 _ yl/2 



M • (6) 



The eigenvalues are in fact bounded and approach zero near the classical singularity fi = 0. 
The fact that the inverse volume is cutoff for small volumes is the principal reason for the 
period of super-inflation for a scalar field in the early universe. We will discuss this more 
later. 

There are two important points we wish to make regarding the inverse volume operator 
given here. First, as noted the formula is obtained using traces of holonomies in the spin 1/2 
representation. The generalization to higher J has been performed and the effect is that 
the cutoff region gets pushed to larger volumes for larger J allowing for phenomenological 



modifications for scales larger than the Planck length (the cutoff for J = 1/2 is deep in 
the Planckian regime where semi-classical equations of motion are no longer expected to 
be valid). Most of the phenomenological investigations are based on the modifications for 
larger values of J. The second point we make is that the holonomy operators (0) involve 
the length parameter /i and an implicit value of //q = 1 has been used to arrive at equation 
(JHJ). As stated in Q| one can use an arbitrary value /xq in the regularization and in that 
reference the Hamiltonian constraint is constructed using a value of v^3/4 determined on 
physical grounds from the smallest allowed area of the full theory of LQG. It is natural that 
the same value should be used both for the Hamiltonian constraint and the inverse volume 
operator since both are regularized using holonomies. We would like to update the formulas 
for the inverse volume operator to account for this ambiguity and in the next section we will 
do so. 

Finally we wish to represent the Hamiltonian constraint operator, classically given in 
equation (j2j). There are two important non-trivialities that arise when attempting to quan- 
tize the constraint. The first is that the connection parameter c is no longer one of the basic 
variables and must be represented using the holonomy operators. The other important non- 
triviality deals with how the sgn(p)y/p term of the constraint gets quantized from the full 
theory. While this term in the isotropic constraint could be simply quantized using the p 
operator jSJ), this approach is not possible in the full theory. There the full constraint is 
given by 3( 

H - l<l\r /VeJ'^F* • ... 



K 



\Q\ 



where N is the lapse, F* b is the curvature of the connection, q is the determinant of the 
three-metric, and the dots indicate the Lorentzian part of the constraint in which we are 
not interested at this point (upon symmetry reduction it is proportional to the first half 
of the constraint). Most important is the presence of the —7=^- term which leads to the 



sgn(p)y/p in the symmetry reduced action, but involves inverse triad operators in the full 
theory. Thus as for the inverse volume this term is represented using a commutator between 
the volume operator and holonomies. We will see that for higher J this terms acquires 
quantum modifications for small volume analogously to the inverse volume operator. 
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The resulting operator is given by 

2% 



GR 



y e ,,fc [hihih7 x h7 v hi 



h k \V 



(7) 



ijk 



where we have explicitly indicated that the trace is performed in the fundamental representa- 
tion. The curvature term of the classical constraint (c 2 ) has been regulated with holonomies 
taken around a closed loop (the h^jh~ l bZ l term). The sgn.(p)y/p term is now represented 



as hi 



h k \V 



in the constraint operator. Several key remarks are to be made. First, the 
holonomy length parameter /xo has been introduced explicitly with the term in the denomi- 
nator and in the holonomy operators. Second, a factor ordering has been chosen in that the 
sgn(p)y/p term has been ordered to the right of the curvature term. This ordering is crucial 
in the singularity removal mechanism as we shall see. However, through this choice in or- 
dering, the constraint is not self-adjoint. Appropriate self-adjoint constraints built from (J7J) 
have been proposed j^ll . 1^ which have been shown to be non-singular. In this paper we will 
consider the non self-adjoint spin J constraint though we note that self-adjoint constraints 
can be constructed in an analogous way which does not effect the results presented here. 

Using this formula we can now determine the action of the constraint operator on the 
basis states Using the formulas for the holonomies in the fundamental representation 
© we get 
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GR 



^2 ~2 

sn cs 



shUcs — csUsh 



whose action on the basis states is 

3 



(V M+W) - V^ ) |/i + 4/i ) - 2\fi) + |/i - 4/x ) 



(8) 



To understand the singularity removal mechanism we need to consider physical wave 
functions which are those annihilated by the constraint operator. Such states can be ex- 
panded using the basis states as \ip) = ^ | yu) . Imposing the constraint equation leads 
to a difference equation for the coefficients if)^ 

3 



ft 



iu+5/io 



V, 



(9) 



where we have assumed that the matter constraint H m acts diagonally on the basis states. 
The key fact about singularity removal is that the wave function coefficient ipo corresponding 
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to the value of the wave functions at the classical singularity decouples from the difference 
equation. This can be seen from the fact that the (V^ — V^_ Mo ) factor that always ac- 
companies ipo on the lhs of the difference equation vanishes. The singularity is then absent 
provided the matter side of the constraint annihilates the \{i — 0) state. This has been 
shown to be true for various forms of matter including a scalar field. 

Let us understand the singularity removal mechanism on a more general level. The key 
observation is that the constraint operator annihilates the \fi = 0) state since (V^^ — V^_ Mo ) 
vanishes identically for /i = in equation (|SJ). Let us assume in general that the constraint 

operator acts by raising and lowering the basis states 

<j 

H\fi)=J2 a l \^ + k ) 

k=—u 

where a and a* are parameters depending on the details of the operator. Assuming this, 
the difference equations for the components ip^ will be 

k= — a 

for all values of //. The statement that the constraint operator annihilates the \\x = 0) state 
implies that = 0. In the difference equation the ipo component is always accompanied 
with the «q term and thus decouples. Had the sga.(p)y/p been ordered to the left of the 
curvature term, the constraint operator would not annihilate the \fi = 0) state and the 
singularity would not decouple. 



III. HOLONOMY REPRESENTATIONS 



We would now like to generalize the calculation of the Hamiltonian constraint operator 
(J7J) and the inverse volume operator to arbitrary irreducible representations labeled by spin 
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J . That this ambiguity is possible in the full theory was first elaborated in 
ticular interest is what effective semi-classical equations of motion can be inferred from the 
constraint operator and what phenomenological differences can arise from the modifications. 
After constructing the arbitrary spin J constraint operator and inverse volume operator, we 
will give an explicit construction of the J = 1 constraint operator and the resulting difference 
equation will be investigated. Following that, we will propose effective classical equations of 
motion and bounds on their validity. 



A. Quantum Constraint Operator and Inverse Volume Operator 

Let us propose the generalization of the formula of (J7J) and show that it is a proper 
representation of the classical Hamiltonian constraint (0). The formula is 

U)H °« = -i^mwjTT) 5 (i>k •M***v* *> ) • (10) 

We can restore the classical constraint by considering the limit as /^o goes to zero. The 
holonomies are given by 

from which we find 

J^^hihjK'hj 1 = 2(fi c) 2 r k + O(fil) 



ij 



h k {h k \V] = -flo{c,V}Tk + 0(f4) =—^KJfio8gn.(p)y/^\T k + 0(j4)- 

Now using the formula ^tr^) = -§ J (J + 1)(2J + 1)8^ we get 

Y,£ lJk Wtr(h i h j hr 1 hJ 1 h k {hi\V})=i4iU(J+l)(2J+l) sgn(p)^\i\c 2 + 0(^) 

ijk 

from which we recover the classical expression given in (J2J). 

As in the previous section the classical expression of the Hamiltonian constraint is pro- 
moted to a quantum operator using the holonomy operators and promoting the Poisson 
bracket to a commutator. The resulting expression is given by 

3i 



,j)S on = ; x>> <** (*m- v*. [*.-■ »] ) 



where we have used the fact that the operator is gauge invariant to arrive at the second line. 
To calculate the action of the operator we need expressions for the matrix elements of the 
holonomies in the spin J representations in order to perform the trace. Fortunately such 
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formulas do exist and a particular form we will use is 

2J-\m+n\ 



(h 



1 1 mn 



T V 

± van / 

s=\m—n\ 
2J-\m+n\ 

s=\m— n\ 

(hi) = e imwi£ 5 

\">3Jmn ° u mn 



y 



cs 



2J-s 



sn 



■\n-m+s 
~ l ) ^2J-s 

— cs sn 



Y 

± IT, 



(13a) 

(13b) 
(13c) 



where the matrix elements m,n G [— J, J] and T mn and Y mn s are constant coefficients given 
by 

T mn = y/(J + m)\(J-m)\(J + n)\(J-n)\ 

1 111 

Y mns = [J + -{m + n - s)]\ [J- -{m + n + s)}\ [-(m-n + s)]! [-(n-m + s)]! (14) 

The index s in the sum increments by two which implies that if (m — n) is an even number, 
the sum over s will comprise even values and vice versa. We can see that the holonomy 
operators are more complicated sums and products of the basic sine and cosine operators 
whose action is given in equations Hence, the holonomies have a well-defined action on 
the basis states 

With the formula for the holonomies (fT3~|) we now can determine the action of the con- 



straint operator. The 
(HU for h 3 . We find 

k 



term is easiest to formulate given the simple expression 



v u rnn u mo v ° 



'tO(J,QC 



V - e tm ^ c V e 



-imuQC 



This operator acts diagonally on the basis states owing to the action of e * m w> c |^ 
|/i — 2m/i ). The eigenvalues are given by 



M = K - y, 



u—2muQ 



(15) 



and so the SU(2) matrix elements of this operator are diagonal. 



Since h\ 



has diagonal matrix elements we need only consider the diagonal el- 



ements ( h\h 2 h 1 l h 2 1 — h 2 hih 2 l h x 



-U-i 



when performing the trace. Using the holonomy 
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formulas ()13all3b|) we get 

2J-\m+n\, 2J-\o+p\, 
j 2J-\n+o\ 2J-\p+rn\ 

( hih,2h 1 1 h 2 1 — h 2 hih 2 1 h 1 1 J = 2_. T mn T no T op T pm N, /_, 

V / mm — ' ' 

n,o,p=-J sl=\m-n\, sl=\o-p\, 

s2=\n— o\ s2=\p—m\ 
\( j\o~n+m~p / ■\n—m+p-o~] 

x(-x) 8 {-l)' M ^ Us 8J - s sn s 

V V V V 

1 mn si 1 no s2 1 op s3 1 pm s4 

8J 

= E z ™ - 8J " S sV 

s=o 

where the parameter S is simply S = si + s2 + s3 + s4 and we have hidden the complexity 
of the rhs of the first line in the parameters which are just constant coefficients. Taken 
together and performing the trace the constraint operator is given by 

J SJ 

( J )ff^ D — — Z S CS 8J_5 sh 5 W - e im Vocy -imu c\ /-i^N 
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1V2J + 1) 

m=-J 5=0 

We can simplify this formula greatly by looking in detail at the coefficients many of 



which vanish. Let us again display the formula 



2J-\m+n\, 2J-\o+p\, 

j 2J— Ira+ol 2J—\p+m\ r , . . , , . n 

■> \( -\o—n+m—p / •\n—m+p—o\ 

3+s4 W l > \ l > J 



Z m — T mn T no T op T pm (*)(!) \- \ \- y- 

7 * mn si I no s2 * op s3 * vmsA 

n,o,p=-J sl=\m-n\, sl=\o-p\, 

s2=\n—o\ s2=\p—m\ 

(17) 

with the restriction that si + s2 + s3 + s4 = S. To simplify this we will use two important 
facts. First, note that when (o — to) + (m — p) in the sum is an even number, the numerator 
vanishes. Second, as we have stated since the s parameters increment by two, their evenness 
is determined by the parameters m, to, o, p, thus for instance if m — n is even then si is even. 
We now show that the coefficients satisfy the following: 



Z s m = for S odd (18a) 
Z° m = ZH = (18b) 
Z s m = -Z s _ m . (18c) 

To show ()18a)) we note that if S = si + s2 + s3 + s4 is odd then this implies that (m — to) + 
(to — o) + (o — p) + (p — m) is also odd. However, this quantity is equal to zero and hence 
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even. Thus S must be even and Z^ = for S odd. (fTFBj) can be shown by calculating 
explicitly. If S = then this implies that all S{ are zero. This only occurs in the sum over 
n, o, p when n = o = p = m. Hence, o — n + m— p = which is even and thus the numerator 
vanishes. A similar calculation holds for . To show (|18c|) one can make the substitution 
in dHD m, n,o,p -> -m, -n, -o, -p. Using the fact that T mn = T_ m „ n and Y mn s = Y_ m _ n s 
(which follow from the definitions ([14)1 ) it is trivial to show that Z^ = —Z_ m . 
Using these properties we simplify the constraint operator to 

o ■ J 4J_1 

( J )j^ R = — V V Z 2S ' cs 8J ~ 2S ' sn 25 ' Jm^y e -imn c (1Q) 

where the term in (fTo^) involving only V vanishes due to ()18c|) . The action of this operator on 
the basis states can be understood as follows. We showed that the term e !mwc l/e" !m '" )C 
acts diagonally on the basis states and in the next section we will show that it approximates 
the sgn(p)y/p term in the classical constraint. The cs and sn terms each raise and lower the 
basis states by discrete steps up to ±8 J. Since the sum is over even powers of cs and sn, H 
will raise and lower the basis states by even values of /xo- Thus in general we can say the 
action on the basis states is 

4J 

^H GR \H)= ^\fi + 2kfi ) 

fc=-4J 

where a 1 ^ are coefficients determined by equation (|19j). The resulting difference equation 
will be of order 8J. Furthermore, the fundamental step size of the difference equation is 
5^ = 2/i . Note that for the J = 1/2 operator some further cancellation occurs and the 
difference equation given in equation © has step size equal to 4/io. While for arbitrary J 
similar cancellations may occur, we will show that for J = 1 the step size is indeed 2//o 
and the difference equation is of order 8J = 8. We can state for certain that the smallest 
possible step size is 2/x . 

With regards to the singularity behavior we have stated that singularity removal occurs if 
the constraint operator annihilates the \fi = 0) state. To show that this occurs for arbitrary 
J we will use the fact that Z^ = —Z^ m . The constraint will thus contain the following term 



m=l 
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whose action on the basis states contributes a J2 m=1 (V r Ai _ 2mA10 — V r M+2mMo ) term. It is easy 
to see that this vanishes for \i = and hence the singularity decouples from the difference 
equation. 

Turning now to the inverse volume operator, the formula for arbitrary J is JscJ 

4z 



v- 



^ (j ; 1)(2J+1) e4* i .^i 



where again we have used gauge invariance to arrive at the second line. Using the formula 



term 



for the holonomy (j!3c|) we find through similar calculations to those of the /13 
of the constraint operator, the inverse volume operator acts diagonally on the basis states 
with eigenvalues denoted by dj 

12 ^ i ^ 

dM = L> j(j + i)( 2 j + i) ]C m V ^°\ (20) 

and can be approximated for large J by the function 

dj(a) « - 3 D\a 2 /al) 
a 6 

D(q) = 2g 1 / 4 [^(g + l) 7 / 4 + sgn(g-l)|g-l| 7 / 4 ]-^^ + l) 11/4 +k-ir i/4 ] (21) 

where a* = yl ~^^ L l p is the characteristic scale factor below which the quantum corrections 
are large. The eigenvalues are bounded with the maximum value occurring near a*. For 
a ^> a* the function dj(a) approximates well the classical expression a -3 . Below a*, dj(a) 
behaves polynomially and can be approximated by treating the sum as an integral to get 

12\V 2 
7 j a? 



^ W « ( - ) ( 22 ) 



Our goal was to note explicitly the role played by the regularization length /i . We note 
that the effect is simply to shift a* to larger or smaller volumes. 

We have thus succeeded in formulating the constraint operator (|19J1 and inverse volume 
operator (J2(Jj) for arbitrary J. The key results for the Hamiltonian constraint are that the 
difference equation is higher order for larger J, yet the step size remains on the same order. 
The singularity removal mechanism remains valid for arbitrary J. For the inverse volume 
operator, the eigenvalues dj are bounded near a critical scale factor a* whose value scales 
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both with /io and J. We do not explore any further the general properties of the quantum 
constraint for arbitrary J as we do not have a simple formula for the coefficients and 
the resulting difference equation would be quite complicated and none the enlightening. 



B. Quantum Theory for J = 1 

We now wish to calculate an explicit example of the quantum operator given in equation 
(fTTjj) . For simplicity we choose the J = 1 representation. We will derive the difference 
equation and show that under appropriate conditions the solutions approximately satisfy 
the Wheeler-DeWitt differential equation. We will also examine the local stability of the 
difference equation to determine the behavior of spurious solutions. 

After a little labor the J = 1 operator is given by 



(1)tj 12 r^6^2 ^4^4 -^2^61 

H-gr — IO q o cs sn — 2cs sn — cs sn 
whence the action on the basis states is 



{1) H GR \fi) = - 3 „ Sl (n) { |M + 8/io) — 4|/x + 6/x ) - % + 4// ) 



32K7 2 yUQ 



+ 4|/i + 2fio) + 6\n) + - 2fio) 
- - 4/i ) - 4|/i - 6/i ) + |/i - 8/i > 



where we have defined 



7^o 

and it can easily be shown that for large volume si(p) ~ -^/p. The resulting difference 
equation is 

QO 3 2 2 < + 8/i ) VV+8 W - 4s i(/^ + 6/io) </V+6 W - 4s i(/ i + 4 )"o) </V+4 M 
32k7 2 /x 

+ 4si (// + 2// ) ^m+2mo + 6s i (A*) ^ + 4si (/i - 2/i ) VV-2w> ( 23 ) 
- 4si(/i - 4/i ) VV-4mo _ 4s i(/ i _ 6 A*o) ^m-6mo + - 8/i ) VV-sw | = -H m {p) ^ • 

As predicted in the general case, the difference equation is of order 8 J = 8 with step size 
8fi = 2//q- 
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To make contact with the Wheeler-DeWitt equation we will make the following assump- 
tions. First is that we are in the large volume limit where \i ^> fiQ. Second is that the wave 
function coefficients ip^ do not vary sufficiently fast on the order of 8/i = 2//o such that we 
can approximate the discrete coefficients with a continuous function ip(fi). Whether or not 
this assumption holds true for large volume is model dependent and an example where this 
continuum approximation fails is the isotropic flat model with a positive cosmological con- 
stant. The consequences for that model are discussed in depth in [24J. Yet, let us assume for 
the moment that we can make these assumptions, hence we can Taylor expand the function 
if}(fj, + 8) ~ if>(fJ>) + + \lT^ 8 2 . Under these assumptions if>(fi) will approximately satisfy 
a differential equation. 

The Wheeler-DeWitt equation is derived from the classical action (J2J) be performing the 
usual Schrodinger quantization where the operator c is quantized to — l^firyj^. The resulting 
differential equation is 

kH 2 d 2 r i 

— [VpHp)\ + Hm^ip) = . (24) 

Turning to the difference equation let us define t M = si(/i)-?/> M and the difference equation 
simplifies to 



— 4t M _ 4At0 — 4t M _6 Mo + t^spo J = —H m (fi) ifjp . 
Using p = and Taylor expanding t{p) we find 

32 2 2 4 <9 2 t(p) 

+ 6t M + 4^_ 2mo - 4t /1 _4 At0 - 4i M _ 6At0 + t M _8 M0 « — g - ^ ^ p Qp2 ■ 

For large volume i(p) = si(p) ^(p) ~ ■Jp'i'ip) an d plugging this in we find that the difference 
equation is approximated by the Wheeler-DeWitt equation. 

Next we would like to examine the local stability of the difference equation. Since the 
difference equation is of higher order than the second order Wheeler-DeWitt equation, we 
must determine the behavior of the spurious solutions. As determined in j^j], if the higher 
order difference equation admits solutions with amplitudes that grow (locally), then the 
difference equation is not locally stable and this might call into the question the validity 
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of the quantization as any semi-classical solutions would quickly become dominated by the 
expanding spurious solutions. 

Let us be more precise. We will consider the large volume limit in a regime where the 
matter contribution is small. We assume that the volume is large enough such that the 
variations of and H m (n) are small. We can thus approximate the difference equation by 
one with constant coefficients 

- 4^_ 4M - 4^_ 6M + ^_ 8M = 32K7 3 g f ( ^ m(//) ^ = P% (25) 

where we assume that on the order of 8ji = 8/i that P = 32K7 3 ^°^ m( ' At ' > is constant. The 
difference equation with constant coefficients can be solved exactly by assuming solutions 
of the form tp^ = Az^^ where z E C. The condition of local stability states that we 
should consider solutions of the homogeneous equation (with P = 0) and that the difference 
equation is locally stable if all solutions have norm equal to one, that is \z\ 2 = 1 [^J. This will 
guarantee that spurious solutions which don't approximate a semi-classical Wheeler-DeWitt 
solution will not come to dominate. 

Plugging in our ansatz ip^ = A z^^° to the homogeneous difference equation (|25J) we find 
the following condition on z 

z * - Az 6 - Az A + Az 2 + 6 + Az' 2 - Az~ 4 - Az~ 6 + z~ 8 = . 

A numerical calculation of the roots shows that there are solutions with norms greater than 
one, thus the J = 1 difference equation is not locally stable. The result of this behavior is 
plotted in figure $1} for a massless scalar field with constant momentum where H m {jj) = 
|dj(//)P|. The figure shows a solution to the difference equation (|2~H|) and a solution to 
the Wheeler-DeWitt equation ()24|) for comparison. Initial conditions are specified on ip^ to 
match the Wheeler-DeWitt solution and the solution is evolved forward using the difference 
equation. It is clear that the solution of the difference equation follows the Wheeler-DeWitt 
solution briefly, after which the spurious non semi-classical solutions quickly dominate the 
wave function. 

The presence of the ill behaved spurious solutions does not by itself represent a problem 
with the quantization. To truly determine the physical consequence requires a detailed 
understanding of the physical inner product. With a physical inner product, one could 
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FIG. 1: Log plot of the difference equation solution (dots) compared with Wheeler-DeWitt solution 
(solid line) for a massless scalar field at large volume. The presence of spurious solutions to the 
higher order difference equation eventually dominate to solution. 

identify which solutions are physical based on the notion of probability, that is unphysical 
solutions would have either vanishing or infinite physical norm and would be modded out 
of the physical Hilbert space. An exampled where this happens is in the symmetry reduced 
Plebanski model of |2J]. There, extra solutions have zero physical norm and the resulting 
physical Hilbert space is one dimensional for matter in the form of a cosmological constant. 
If the spurious solutions here do indeed have zero physical norm then the quantization 
presented would have the correct semi-classical limit. If not this would indicate the need 
for different quantizations or even the possibility that the higher spin quantizations do 
not represent consistent quantizations of isotropic LQC This might indicate a preference 
of the quantum theory toward the use of the fundamental representation for defining the 
Hamiltonian constraint. There are possible indications that this behavior also occurs in 
the full theory where a large set of spurious solutions exist for the higher J gravitational 
constraint operator j^J. 
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C. Effective Classical Theory 



We now turn to the question of what sort of semi-classical equations of motion can be 
inferred from the constraint operator given in equation (|19j) . While in general the effective 
equations of motion should be derived from the wave function solutions themselves (for in- 
stance by finding semi-classical states), here we will take the simplest route. Our motivation 
comes from the path integral quantization of LQC where a discretized path integral involves 
integration over a classical action (for details on the path integral derivation see A 
detailed derivation of the path integral shows that the effective classical Hamiltonian con- 
straint is given by the p-q symbol H e ff = ^ C J?J,^ . Since the Hamiltonian constraint of 
LQC consists of a self-adjoint part (the hihjh~ x h^ x term) of which |c) is an eigenstate, and 
a part [h k {h^ 1 , V}) of which \p) is an eigenstate we need only consider the eigenvalues 
of the two operators to get an effective classical constraint written in terms of a classical 
c and p. We note that a more detailed consideration involving semi-classical states could 
lead to additional corrections. Equations of motion can then be derived from the effective 
Hamiltonian through the Hamiltonian equations p = {p,H e ff} and c = {c, H e ff}. We are 
interested in determining what modifications arise for instance to the Friedmann equations 
from the quantum effects. 

We now show that two kinds of corrections arise in the effective constraint. The first 
occurs for small volumes and arises from the hk {h^ , V} term which we will show behaves 



as a modification to the classical sgn(p) yjpf term in the constraint. The second modification 

is due to the curvature term h;h~h~ x h~ x and amounts to a modification of the classical c 2 
term. This modification is evident when the connection c is large or equivalently (for flat 
models) when the extrinsic curvature is large. We will show that for large volume and small 
extrinsic curvature, the correct classical equations are recovered. 

Let us now calculate the corrections explicitly. From equation (JT5J) we get the formula 
for the h 3 {h 3 1 , V} term as 





mn 




mn 



Next the curvature term hihjh { hj can be calculated by expanding the holonomies hi 
e vocn ag done for equation pip. If taken to order (/xqc) 4 we find that 



h x h 2 h- x x h~ 2 x - h 2 h 1 hz 1 h- 1 1 = 2(fi c) 2 T 3 + (/ioc) 3 (7i - r 2 ) - -(/i c)V 3 + C((/i c) 5 ) . 
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Now using (r 3 ) mn = im5 mn and the fact that (ri) mn and (r 2 ) mn are off-diagonal matrices we 
find after performing the trace 

^ fj,— 2mno\ 



,3^5, 



-,sj c*-UA+0(^) (26) 



K7 2 \ 3 / V K7 2 

and this is the formula we will use as our effective constraint. Notice the similarity to the 
classical constraint in (j2J) except for two modifications. The first is the presence of the 
function sj which we have defined to be 

6 J 

SJ = ~ 7 Z>J(J + 1)(2J+1) m 5 j mV2mM (2?) 

which we now show approximates sgn(p)y]p|. This can be seen by noting that for large 
volume V^ 2mm « - m lf i l 2 p /2 sgn(p) y/\p\ + O(^). Thus we get 

~ sgn(p)v^N + °(^^ 

where we have used J^m=-J m2 = ^(^ + tyftJ + 1) to arrive at the second line. As for 
the formula for dj{fx) we can approximate sj(/j,) for large values of J by treating the sum 
as an integral. The resulting formula is 



it 2^piir + 2J/io)5/2 + sgn(/i ~ 2J/io)l/i ~ 2J/io|5/2] 

^[( /U + 2J/io) 7/2 -|/x-2J/i | 7 / 2 ]} 



35 

which implies that written in terms of the scale factor a 
sj(a) = aS{a 2 /al) 

%) = ^{^[(s + l) 5/2 + sgn(g-^ 



where as for dj(a) we have the same critical scale factor a* = y -^f^-lp below which quantum 
corrections occur. The function S(q) for q > 1 is approximately equal to one which implies 
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that for a > a*, sj(a) m a — \/]p\- For small volume S(q) behaves as \y/q hence sj(a) 
behaves quadratically with a 

sj(a)^- — a < a* . (29) 
5 a* 

A plot of the function sj(a) compared with ^fp = a is shown in figure (J2J). It is clear that sj 
behaves quadratically for small a and changes behavior near a* after which it approximates 
the classical expression ^fp. 




a 



FIG. 2: LQC modified function sj(a) (|27j) compared with the classical expression for the function 
sgn.(p)y/\p\ = a. The modified sj(a) function behaves quadratically for a < a* = ^5%. The 
values of a* in the graph are 1.290, 2.89, 4.08 for J = 10, 50, 100 respectively. 

Thus we have determined the effective Hamiltonian constraint given in equation ()26|) 
which is valid in the regime where the connection component c is small or more precisely 
when /zo c ^ 1- The function sj(/i) in turn approximates the classical sgn(p)^/\p\ term for 
large volume when \i 3> 2J/z or equivalently when a ^> a*. 

Let us now make contact with the standard formulation of isotropic cosmology written 
in terms of the scale factor a. We assume that the orientation of the triad is positive 
(p, \x > 0) so we can drop the absolute value and sign terms in the formulas. According to 
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equations pal |3J) we have the relation between p and the scale factor as a 2 = \pl 2 v - From 
the Hamiltonian equations and the vanishing of the constraint we can derive the Friedmann 
equation. We first have 

• r u \ 1 dH 'ff 
= l(2c-~p 2 c 3 )sj. 

Using p = a 2 we can use this equation to get the Friedmann equation in terms of the 
connection as 

«V S J° 2 ( ^ 4 2 2 4 4 4 

From this it is clear that the criteria poc <C 1 corresponds to a criteria on the scale factor 
velocity being small po'J— o "C 1. To get this equation in standard form we need to solve 
for the connection in terms of the matter Hamiltonian. We can get this from the vanishing 
of the full Hamiltonian constraint H = H e ff + H m = 0. This gives 

~ 3s j m+ 27s 2 m 

where we have only kept the terms to second order in H m . Putting these two together we 
get the modified Friedmann equation 

a\ 2 K Sj K 2 7Vo tj2 



a) 3 a 4 9a 4 

= ^ S(a) p m -^K 2 j 2 p 2 Q a 2 p 2 m (30) 

where the matter density is simply p m = a~ 3 H m and as before sj(a) = a S(a). The modified 
Friedmann equation remains valid provided p^—a 1. This in turn places a bound on 



the matter density p m 3 f ^2 below which the effective equations remain valid. From the 
modified Friedmann equation it is clear that in the large volume regime (where S(a) ~ 1) 
and small curvature (where we neglect the p 2 m term) we recover the standard Friedmann 
equation (f) 2 = fp m . 

Let us reiterate the main results of the section. The modified effective Hamiltonian 
constraint is given in equation ()26|) and we rewrite it here 

Heff = -^- 2 sj (c 2 - \p 2 A (31) 



«7 2 \ 3 
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where the function s,j is given in equation (}2*Tj) and approximates sgn(p)^^/\p\ for large 
volume. The modified Friedmann dynamics given in equation ()30|) was calculated through 
the Hamiltonian equations derived from the effective constraint. The effective constraint 
is valid when /i c <C 1 which in turn places a bound of validity determined by both the 
extrinsic curvature and on the matter density. 

IV. DISCUSSION 

In this paper we have sought to clear up certain issues in LQC, namely the complete 
consequences of quantizing the Hamiltonian constraint using higher spin representations for 
the holonomies. In addition we have fixed the constants necessary to make contact on a 
consistent basis with the standard action of GR. We have found that the use of higher spin 
holonomies to regulate the gravitational part of the constraint operator leads to modifications 
of the det(|g|) term. The modifications are qualitatively similar to those of the inverse scale 
factor which lead to the dj(a) function. This is not unexpected since we are quantizing the 
constraint operator to take into account factors involving the inverse triad (the l/det(\q\) 
term). This modification, leading to the sj(a) function, is another important consequence 
of a key mantra of LQC: to remain as close to the full theory of loop quantum gravity as 
possible. We could have exploited the symmetry of the isotropic model and quantized the 
sgn(p)yjp| term avoiding any inverse triad operators. Because this is not possible in the 
full theory it is not done in LQC and the non-perturbative corrections are imported from 
the full theory. 

We have been guided by certain simplicity considerations. Previous studies have used the 
J = 1/2 gravitational part of the constraint (and the resulting difference equation) while 
freely choosing an arbitrary representation to define the inverse volume operator in the 
matter part of the constraint. While there is certain freedom to specify the representations 
differently and quantize the two parts in a different manner, for instance to define a Jq for 
the gravitational part and a J m for the matter part, the simplest choice is to use the same 
representation. This ambiguity also exists in the parameter /i appearing in the matter and 
gravitational part. 

Thus it is important to reexamine the phenomenological consequences in light of the 
modifications to the Friedmann equation ()3U|) that arise for arbitrary values of J . Let us 
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first concentrate on the modifications arising from sj(a). In the instance of a massless scalar 
field we have for the Friedmann equation 



where we have used H m = 7}dj(a)P? with P^ being constant. To understand the behavior 
of the scalar field let us de-parameterize the equations of motion to remove the notion of 
time. This is in accord with notions of quantum gravity where coordinate time does not 
have any physical meaning. Thus we are interested in how the scalar field evolves with the 
scale factor as opposed to coordinate time. From the Hamiltonian equations we find that 
= dj{a)P ( j > and we de-parameterize as d(f)/da = <p/a to get 



which implies that the field is pushed to higher values in the region below a* since the S(a) 
appearing in the denominator is suppressed in that region. The addition of the s j corrections 
can therefore improve the viability of LQC to push the inflaton up its potential to set the 
initial conditions for slow-roll inflation. 

The consequences of the modifications to the Friedmann equation quadratic in matter 
have not been explored. Important though is that we have given a precise bound a^—a <C 1 

S j 

on the validity of the first order effective constraint. This is in contrast with the criteria for 
validity cited in other phenomenological investigations a/o> y/^lp 0] based on the notion 
that the Hubble length should not be smaller than the fundamental length of discreteness 
■J^l v . The criteria given here can be understood heuristically as follows. The discreteness 
of the difference equation plays a role precisely where the Wheeler-DeWitt approximation 
f!24j) breaks down, that is when the wave function changes significantly on the order of the 
fundamental step size Sfi = 2/i . If we consider a Wheeler-DeWitt solution that is locally 
oscillatory ip(^) = exp(i2ix fx j 'A) where A is the wavelength, then the discreteness becomes 
important when A ~ 4/i . Given that in the Wheeler-DeWitt equation c = ^-^-d/dp we 
find that the connection c is related to A as c = 47r/A. Using this we recover the criteria 
that the higher order corrections occur when /ioc ~ 1. We can therefore understand that 
these corrections are a direct result of discreteness effects of the difference equation. 

The term quadratic in the matter density is negative definite in the Friedmann equation 
and hence can act as an effective compactness even in the non-compact model considered 
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in this paper. Whether or not this would indeed imply a re-collapse depends on the in- 
clusion of the full corrections to arbitrary order of [IqC. Furthermore, classical descriptions 
in this regime are not expected to be valid since, as stated in the previous paragraph, this 
corresponds to the regime where the discreteness is dominant. A model where this occurs is 
deSitter space which is discussed in detail in j^J. In that model classically the extrinsic cur- 
vature grows with time and the wave functions become non semi-classical around a critical 

( \ 3//2 

volume given by V c = ( ^2^2 J which is precisely when fi c = ir. Above the critical vol- 
ume, the wave function decays rapidly indicating a classically forbidden region and thus an 
effective compactness. Thus this correction can have drastic consequences for the evolution 
of the universe. For ordinary matter such as dust or radiation, the matter density drops off 
sufficiently fast for these corrections not to be relevant at large volumes. Yet the fact that 
observations suggest that the universe might be in an asymptotically deSitter space call into 
question the validity of LQC for large volumes and a more detailed derivation of LQC from 
the full theory would likely solve this problem. 

We are in the process of examining numerically the full phenomenological consequences 
of the modified Friedmann equation (jHOj) for more complicated forms of matter. The fact 
that the sj corrections help push the inflaton up its potential is a benefit to this scenario 
and would help enlarge the parameter space that leads to successful inflation. 

The key parameter involved in the quantum corrections is the regulating length /a . A 
deep understand of the origin of this parameter is essential before exact phenomenological 
predictions can be made. In the full theory one takes the limit as /io goes to zero and 
a well-defined operator remains which due to diffeomorphism invariance is independent of 
the regulator 19(. In LQC the corrections for sj, dj as well as for large extrinsic curva- 
ture depend explicitly on the value of /io . If we were to take the limit of small /io, the 
phenomenological consequences of the quantum corrections would be pushed into the deep 
Planckian regime. This might wash out any possible observational signals from LQC. Yet, 
a smaller value of /io would also push to larger volumes the validity of LQC for the deSitter 
example discussed above. In lj| it is argued heuristically that /io is fixed to a value of \/3/4 
in the Hamiltonian constraint based on the smallest allowed area of the full theory. A more 
thorough understand of the parameter /x and its derivation from the full theory is required 
to settle this issue. 

A key issue with which we have not dealt is the need to derive the effective equations of 
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motion from solutions of the difference equation itself. Several proposals have been put forth 



to date on this matter. In one 



2l\ , kinematical semi-classical states are constructed and the 



effective constraint is calculated from the expectation value of the constraint operator. The 
validity of this is not entirely clear since the semi-classical states are not physical states in 
that they are not annihilated by the constraint operator. The corrections to the constraint do 
include a /igC 4 term as we showed in this paper (in that paper the corrections are calculated 
for large volume for the J = 1/2 constraint so would not include sj corrections). The exact 
factors for these corrections do not agree however. In addition corrections arise in that 
model due to the spread of the wave function, something of which we have not considered 
here in deriving the effective equations of motion. 

In another technique 2^, WKB type solutions of the difference equation are calculated 
and from that an effective constraint is extracted. In this paper the modified constraint 
includes a —-^r^ [ s i/2(/W + 4/io) + Sx/ 2 (/i — 4/i )] c 2 term which to leading order is equivalent 
to the effective constraint of equation (12611 for J = 1/2. More detailed calculations in 

show the higher order corrections in /xqc appearing in our effective 
constraint 1)26)1 . In that paper the effective constraint for large volume is given by H e ff = 
— K7 2^ y/p sin 2 (/xoc) (where we have adjusted the multiplicative factors to agree with ours) 
which is precisely the effective constraint presented in this paper had we calculated all the 
higher order corrections for the J = 1/2 constraint (the c 2 — \^c i term can be seen as the 
small /ioc expansion of sin 2 (/xoc)//Xo). In addition the WKB effective constraint contains a 
potential term which contributes for small volume. Again the interpretation of dynamics 
in the WKB setting is not entirely clear especially in the setting of the discrete difference 
equation of LQC In addition the approximation used in 28] relies on correlating wave 
function solutions of the difference equation on scales smaller than the fundamental step 
size Sn = 4/i . It can be shown in LQC that quantum interference only occurs for quantum 
states defined at values of \i differing by an integer times the fundamental step size 3|. 
Thus, quantum mechanically LQC does not correlate the quantum wave function at a given 
volume ji and a neighboring value /x + e where e < 5ji. 

A further method to determine semi-classical dynamics involves introducing a coordinate 



221. In that 



time parameter with which to consider evolution of Gaussian wave packets 
paper the wave packets are shown to follow the classical trajectory even to small volumes. 
The relevance of the solutions to physical solutions (annihilated by the constraint operator 
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and independent of any time parameterization) has yet to be determined. 

Each of the techniques to derive effective equations of motion are performed with models 
with zero degrees of freedom, and thus none consider physical semi-classical states. In 
quantum gravity where coordinate time is not a physical degree of freedom, one needs to 
include higher degrees of freedom and choose one of them to play the role of a clock. The 
model of a massless scalar field provides a testing ground for these ideas. In this model, 
the volume is a monotonically increasing function (classically) and can play the role of a 
clock. A semi-classical state peaked around some value of the scalar field at a given volume 
can be evolved forward in "time" with the difference equation and the resulting trajectory 
compared with the effective equations of motion. This would have the benefit that the wave 
functions considered would be physical states annihilated by the constraint operator. Early 
results for the spin 1/2 difference equation indicate that the first order corrections for the 
effective equations match well with the results from the difference equations to volumes near 
the Planck scale. While this issue would be trickier to study for the higher order spin J 
difference equation, the agreement for the spin 1/2 constraint is indicative of the validity of 
the effective equations. 
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